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We propose a method which we call "Isotropic Entanglement" (IE), that predicts the eigenvalue 
distribution of quantum many body (spin) systems (QMBS) with generic interactions. We inter- 
polate between two known approximations by matching fourth moments. Though, such problems 
can be QMA-complete, our examples show that IE provides an accurate picture of the spectra well 
beyond what one expects from the first four moments alone. We further show that the interpolation 
is universal, i.e., independent of the choice of local terms. 

PACS numbers: TBF 



A HIGHLY ACCURATE MATCH (BETTER 
THAN FOUR MOMENTS). 

We propose a method to compute the "density of 
states" (DOS) or "eigenvalue density" of quantum spin 
systems with generic local interactions More gen- 
erally one wishes to compute the DOS of the sum of 
non-commuting random matrices from their, individually 
known, DOS's. 

We begin with an example in Figure [1] where we com- 
pare exact diagonalization against two approximations 
that we considered early in our work and our method 

• Dashed grey curve: classical approximation. No- 
tice that it overshoots to the right. 

• Solid grey curve: isotropic approximation (or iso). 
Notice that it overshoots to the left. 

• Solid black curve: isotropic entanglement (IE). 

• Dots: exact diagonalization of the quantum prob- 
lem given in Eq. [TJ 

The classical approximation ignores eigenvector struc- 
ture by summing random eigenvalues uniformly from 
non-commuting matrices. The curve is the convolution 
of the probability densities of the eigenvalues of each ma- 
trix. 



N=5, d=2, r=4, trials=500000 
Isotropic Approximation 



Isotropic Entanglement (IE) 




Figure 1: The exact diagonalization in dots and IE com- 
pared to the two approximations. The title parameters are 
explained in the section on numerical results. 



The isotropic approximation assumes that the eigen- 
vectors are in "general position"; that is, we add the two 
matrices with correct eigenvalue densities but choose the 
eigenvectors from Haar (uniform) measure. As the ma- 
trix size goes to infinity, the resulting distribution has 
been popularized recently as the free convolution of the 
individual distributions [11, or 0, for a nice exposition]. 

The exact diagonalization given by red dots, the 
dashed and solid grey curves have exactly the same first 
three moments, but differing fourth moments. 

Isotropic Entanglement (IE) is a linear combination of 
the two approximations obtained by matching the fourth 
moments. We show 1) the fit is better than what might 
be expected by the first four moments alone, 2) the com- 
bination is always convex for the problems of interest, 
given by a < p < 1, and 3) this convex combination is 
universal depending on the parameter counts of the prob- 
lem but not the eigenvalue densities of the local terms. 

Parameter counts: exponential, polynomial and 
zero. — Because of the locality of generic interactions, the 
complete set of eigenstates has parameter count equal to 
a polynomial in the number of spins, though the dimen- 
sionality is exponential. The classical and isotropic ap- 
proximations have zero and exponentially many random 
parameters respectively. This suggests that the problem 
of interest somehow lies between the two approximations. 

PHYSICAL IMPLICATIONS AND PROBLEM 
FORMULATION. 

QMBS spectra have been elusive for two reasons: 1. 
The terms that represent the interactions are generally 
non-commuting. This is particularly pronounced for 
systems with random interactions (e.g., quantum spin 
glasses i, p. 320][4, ?!). 2. Standard numerical diag- 
onalization is limited by memory and computer speed 
because of exponential growth. Energy eigenvalue dis- 
tributions are needed for calculating the partition func- 
tion Z . Exact calculation of the spectrum of interacting 
QMBS has been shown to be difficult [6]. 

Though, much progress has been made in understand- 
ing the ground states of interacting QMBS [3-[l3|> eigen- 
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value distributions are less studied. An accurate descrip- 
tion of tails of distributions are desirable for condensed 
matter physics. IE provides a direct method for obtaining 
eigenvalue distributions of quantum spin systems with 
generic local interactions (i.e., quantum spin systems) 
and does remarkably well at approximating the tails. 

Though we are not restricted to one dimensional 
chains, for sake of concreteness, we investigate N inter- 
acting d-dimensional quantum spins (qudits) on a line 
with generic interactions. The Hamiltonian is 



Af-l 

1=1 



i?L---,i + L-l <X)IdN-!-(L-l), (1) 



where the local terms Hi ... are finite x 

random matrices. From now on we take the case of 
nearest neighbors interactions, L — 2, unless otherwise 
specified. 

The eigenvalue distribution of any commuting subset 
of H such as the terms with I odd (the "odds") or I even 
(the "evens") can be obtained using local diagonalization. 
However, the difficulty in approximating the full spec- 
trum oi H = H^^^ + ^even ^® summing the odds and 
the evens because of their overlap at every site. 

The intuition behind IE is that terms with an overlap, 
such as and iJ;+i.j+2, introduce randomness and 

mixing through sharing of a site. Namely, the process of 
entanglement generation introduces an isotropicity be- 
tween the eigenvectors of evens and odds that can be 
harnessed to capture the spectrum. 

ISOTROPIC ENTANGLEMENT. 



The eigenvalue distribution of Eq. |3] is what we are 
after. 

To establish the plausibility of the classical and 
isotropic approximations, we consider the moments. In 
general we have the k*"^ moments 

ml'" ^-^ETi- {A + Q-^BQ)'' and (5) 

The first three moments are usually encoded as the 
mean, variance, and skewness; the fourth moment is 
encoded by the kurtosis which we denote 72 . 

(The Matching Three Moments Theorem) The 

first three moments of the classical and isotropic approx- 
imations exactly match those of the quantum problem 

Turning to the fourth moment, we propose to match 
the kurtosis 7I with a linear combination of the classical 
(7I) and isotropic (72"°) kurtoses: 



7|-P72 + (1-P)7r 



P=^-l^- (6) 

72 - 72 



In terms of probability measures, IE provides 

« dv'^ = pdv" + (1 - p)dv''°, (7) 

where dv denotes a probability measure. So long as 
< p < 1, dv^^ is a probability measure whose first 
four moments match the theoretical measure dv'' . 



Let A and B denote diagonal matrices containing the 
eigenvalues, in random order, of the odds and evens re- 
spectively. In a basis where (at least) the odds are di- 
agonal, our approximations to the Hamiltonian in Eq. [T] 
become 



Classical: H" = A + B 
Isotropic: ^ A + Q^^BQ, 



(2) 



where Q is a /?— Haar measure orthogonal matrix of size 
with (3 = 1 corresponding to real orthogonals. f3 — 2 
unitaries, and /3 — A quaternions (see for formal 
treatment of a general /3). 

The Hamiltonian given by Eq. [1] in a basis where the 
odds are diagonal is 



Quantum: H = A + Q-^BQq. (3) 

The Qq has the formidable form Qq ~ {Qa) ^ Qb, 
where (for odd number of sites N) 

Qa =Qi «) Q3 «> • • • «) Qn-2 ® Id 

QB^h<»Q2<»Q4<»---®QN-i- (4) 



In the expansion of the fourth moments, by a theorem 
we call The Departure Theorem[l|, the numerator and 
the denominator in Eq. IH] respectively become. 



7i - ir 



72 - 7?° 



2 E {Tr [{AQ-^BQqY - (AQ-^BQ)^] } 



d^ 



2 E{Tr[(^B)2-(AQ-iBQ)']} 



d^ 



(8) 



Evaluation of p in Eq. [S] reduces to the evaluation of 
the right hand sides of Eq. [51 It is natural to ask, does a 
< p < 1 always exist such that the convex combination 
Eq. [5] can be formed? How would quantum spectra look 
in the thermodynamical limit (N — > oo)7 The Slider 
Theorem provides the answers. 

(The Slider Theorem) The quantum kurtosis lies in 
between the classical and the isotropic kurtoses, 72"° < 
72 < 72. Therefore there exists a < p < 1 such that 
72 ~ P72 + (1 ^ p) 72''°- Furthermore, lim^v^oo P ~ ^■ 

Through our numerical investigations, we observed 
that p did not change by the choice of local terms; it did 
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not dependent on the covariance matrix. This lead us to 
prove that p was universal with respect to the choice of 
generic local terms and allowed us to derive an analytical 
formula (Figure [21 where we take 13 ~ 1). 

(Universality Corollary) p {N, d, /?) , namely, it 
is independent of the distribution of the local terms (see 



N=3, d=5, trials=50000 




Figure 2: An example where /3 — 1: the quantum problem for 
all d lies in between the iso {p — 0) and the classical {p — 1). 
Entanglement shows itself starting at the fourth 
moment; further, in the expansion of the fourth 
moments, only the terms that involve a pair of local 
terms sharinn a site differ 

N=1 1 , dS; r=4, trials=2800 



— p = 

■ --p = 1 

• Exact Diagonalization 
■ E. Theory: p = 0.73 




Figure 3: A*' = 11 with full rank Wishart matrices as local 
terms. 



N=6, d=4. L=5, r=1024, trials=14 



• Exact Diagonalization 
— i.E. Theory: p = 
-'-'Standard Fit 




2000 3000 
Eigenvalue 



Figure 4: Approximating the quantum spectrum by H^^ — 
Y^i^iQi Hi,-- ,1+aQi, where the interaction is taken to be 
5— local Wishart matrices. 



NUMERICAL RESULTS. 

Here we compare our theory against exact diagonaliza- 
tion for various number of sites N . local ranks r, and site 







---p = 

p = 1 

• Exact Diagonalization 
— IE Theory: p =0.046 
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Eigenvalue 

Figure 5: Local terms have a random binomial distribution. 

dimensionality d. As our first two examples (Figures [T] 
and[3]) we take the local terms to be independent Wishart 
matrices: Each = A^A, where A is a. d^ x d^ ma- 

trix with independent and identically distributed Gaus- 
sian random entries. The higher moments of the Wishart 
matrices were obtained using MOPs 



i5i. 



One would expect that L > 2 would be closer to the 
isotropic approximation than when the interactions are 
nearest neighbors because of the larger number of ran- 
dom parameters in Eq. [TJ When the number of random 
parameters of the problem {N — L + 1) d^ becomes com- 
parable to d^ , we indeed find that we can approximate 
the spectrum with a high accuracy by taking the sum- 
mands to be all isotropic (Figured]). 

Lastly take the local terms to have Haar eigenvectors 
but with Bernoulli eigenvalues: = QfMQi, where 

A; is a diagonal matrix of random eigenvalues ±1 (Figure 
[SJ . Here classical treatment leads to a binomial distribu- 
tion. As expected p — 1 in Figure \5\ has three atoms 
at —2,0,2 corresponding to the randomized sum of the 
eigenvalues from the two local terms. The exact diago- 
nalization, however, shows that the quantum chain has 
a much richer structure closer to iso; i.e, p = 0. This is 
captured quite well by IE with p = 0.046. 

Most distributions built solely from the first four mo- 
ments would give smooth curves. Roughly speaking, the 
mean indicates the center of the distribution, variance 
its width, skewness its bending away from the center and 
kurtosis how tall and skinny versus how short and fat 
the distribution is. It is hard to imagine that the kinks, 
cusps and local extrema of the quantum problem (as seen 
in some of our examples and in particular Figure[5]) could 
be captured by fitting only the first four moments of the 
QMBS Hamiltonian to a known distribution. Remark- 
ably a one parameter (i.e., p) interpolation between the 
isotropic and classical suffices in capturing the richness 
of the spectra of QMBS. Figure [6] summarizes IE. 

ACCURACY BEYOND FOUR MOMENTS. 

We illustrate in Figures [7| and [8] how the IE fit is better 
than expected when matching four moments. We used 
the exact first four moments to approximate the density 
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add random variables = 
Classical convolution of densities 



add random variables = 
Classical convolution of densities 



Several options for adding 
random variables 



Isotropic (or free) convolution = 
dv'^° = dv. ffl,,„dv„ 



Isotropically Entangled = 
(1-p) dy'^^H- p dv"^ 



Classical convolution ^ 
dv^ = dv, 1 dv„ 



Isotropic 
p = 



Quantum convolution = dv'' = dv^ ^q^B 
Isotropically Entangled 



Classical 
p = l 



Figure 6: Summary of Isotropic Entanglement. We have de- 
noted the known classical and isotropic convolutions by * and 
ffljgQ respectively. IE approximates the unknown quantum 
convolutions ffl„. 



4-moment matching mettiods vs. IE: binomial local terms 





• Exact Diagonalization 






— I.E. Theory: p = 0.046 






— Pearson 








— Gram-Charlier 













Eigenvalue 

Figure 7: IE vs Pearson and Gram-Charlier fits for local terms 
with binomial distribution 



using the Pearson fit as implemented in MATLAB and 
also the well-known Gram-Charlier fit In it was 
demonstrated that the statistical mechanics methods for 
obtaining the DOS, when applied to a finite dimensional 
vector space, lead to a Gaussian distribution in the low- 
est order. Further, they showed that successive approxi- 
mations lead naturally to the Gram-Charlier series 



Comparing these against the accuracy of IE leads us to 
view IE as more than a moment matching methodology. 

OUTLOOK. 

Our work supports a very general principle that one 
can obtain an accurate representation of inherently ex- 
ponential problems dealing with QMBS by approximat- 
ing them with far less complexity. This realization is at 
the heart of other recent developments in QMBS research 
such as Matrix Product States [8, 9], and Density Matrix 
Renormalization Group [lo| , where the state (usually the 
ground state of ID chains) can be adequately represented 
by a Matrix Product State (MPS) ansatz whose param- 



4-moment matching methods vs. IE: Wishart local terms 



• Exact Diagonalization 

I.E. Theory: p = 0.42817 

— Pearson 

Gram-Charlier 




Figure 8: IE vs Pearson and Gram-Charlier fits for local terms 
with Wishart distribution 



eters grow linearly with the number of quantum parti- 
cles. Future work includes explicit treatment of fermionic 
systems and numerical exploration of higher dimensional 
systems. 
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